Variational Approach to Hydrogen Atom in Uniform 
Magnetic Field of Arbitrary Strength 



O 

o 
o 

^ : 

a: 



M. Bachmann, H. Kleinert, and A. Pelster 
Institut fur Theoretische Physik, Freie Universitat Berlin, Arnimallee 14, 14195 Berlin 

(February 1, 2008) 

Extending the Feynman-Kleinert variational approach, we calculate the temperature-dependent 
effective classical potential governing the quantum statistics of a hydrogen atom in a uniform mag- 
netic at all temperatures. The zero-temperature limit yields the binding energy of the electron which 
is quite accurate for all magnetic field strengths and exhibits, in particular, the correct logarithmic 
growth at large fields. 



I. INTRODUCTION 
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The quantum statistical and quantum mechanical properties of a hydrogen atom in an external magnetic field 
are not exactly calculable. Perturbative approaches yield good results only for weak uniform fields as discussed in 
detail by Le Guillou and Zinn- Justin , who interpolated with analytic mapping techniques the ground state energy 
between weak- and strong-field. Other approaches are based on recursive procedures in higher-order perturbation 
theory . Zero-temperature properties were also investigated with the help of an operator optimization method in 
a second-quantized variational procedure Q . The behaviour at high uniform fields was inferred from treatments of the 
I , one-dimensional hydrogen atom . Hydrogen in strong magnetic fields is still a problem under investigation, since 
its solution is necessary to understand the properties of white dwarfs and neutron stars, as emphasized in Refs. |^| p| . 
■^j- , A compact and detailed presentation of the bound states and highly accurate numerically values for the energy 
C***- ■ levels is given in Ref. ]l2|| . 

Equations for a first-order variational approach to the ground-state energy of hydrogen in a uniform magnetic field 
based on the Jensen-Peierls inequality were written down a long time ago Jl3[ , but never evaluated. Apparently, they 
merely served as a preparation for attacking the more complicated problem of a polaron in a magnetic field []1 3| — |l 5f| . 

In our approach, we calculate the quantum statistical properties of the system by an extension of variational 
perturbation theory []l6|| . The crucial quantity is the effective classical potential. In the zero-temperature limit, it 
yields the ground state energy. Our calculations in a magnetic field require an extension of the formalism in Ref. |l6| ] 
which derives the effective classical potential from the phase space representation of the partition function. 

Variational perturbation theory has an important advantage over other approaches: The calculation yields a good 
effective classical potential for all temperatures and coupling strengths. The quantum statistical partition function is 
obtained from a simple integral over a Boltzmann-factor involving the effective classical potential. The ground state 
^J^" energy is then obtained from its zero-temperature limit. The asymptotic behaviour in the strong-coupling limit is 
^ ' emerging automatically and does not have to be derived from other sources. 
•<-h , 

II. EFFECTIVE CLASSICAL REPRESENTATIONS FOR THE QUANTUM STATISTICAL PARTITION 

FUNCTION 

A point particle in D dimensions with a potential V(x) and a vector potential A(x) is described by a Hamiltonian 
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#(P> X ) = 7T77 P--A(x) +V(x). (2.1) 



The quantum statistical partition function is given by the euclidean phase space path integral 

Z = j V' D xV D p e - A[p '* ]/n (2.2) 

with an action 

A[p,x}= drHp(r).x(T)+%(r),x(r))], (2.3) 
Jo 

and the path measure 
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V' D xV D p= lim 

JV->cx 



N + l r 

n 

n=l 



(2-Kh) 
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(2.4) 



The parameter /? = l/fc^T denotes the usual inverse thermal energy at temperature T, where fcs is the Boltzmann 
constant. From Z we obtain the free energy of the system: 
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-In Z. 



(2.5) 



In perturbation theory, one treats the external potential V(x) as a small quantity, and expands the partition function 
into powers of V(x). Such a naive expansion is applicable only for extremely weak couplings, and has a vanishing radius 
of convergence. Convergence is achieved by variational perturbation theory |l6| ], which yields good approximations 
for all potential strengths, as we shall see in the sequel. 



A. Effective Classical Potential 



All quantum-mechanical systems studied so far in variational perturbation theory were governed by a Hamiltonian 
of the standard form 



tf(p,x) = ^ + nx). 



(2.6) 



The simple quadratic dependence on the momenta makes the momentum integrals in the path integral fl2.2|) trivial. 
The remaining configuration space representation of the partition function is used to define an effective classical 
potential \4ff(xo), from which quantum mechanical partition function is found by a classically looking integral 



Z = 



d u x 



exp [-/3V c s(xq)] , 



(2.7) 



where Ath = y 2irTi 2 [3/M is the thermal wavelength. The Boltzmann factor plays the role of a local partition function 
Z x ° , which is calculated from the restricted path integral 



-/3V eff (x ) 
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<j)V D xS{x -x(T))e- 



A[x]/h 



with the action 



and the path measure 



A[x] = / gIt 
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-x 2 (r) + y(x(r)) 



lv D x = lim TT ( f 



d D x n 



h 2 (3/M(N + 1)] D ' 2 



The special treatment of the temporal average of the Fourier path 

r hf3 



1 f np 

x Q =x(r) = -y o drx(r) 



(2.8) 



(2.9) 



(2-10) 



(2.11) 



is essential for the quality of the results. It subtracts from the harmonic fluctuation width (x 2 ) cl the classical 
divergence proportional to T = l/ksP of the Dulong-Petit law [|l6|jl§| ]. Such diverging fluctuations cannot be treated 
perturbatively, and require the final integration in expression ( p77| to be done numerically. 

For the Coulomb potential V(x) = — e 2 /Aite^ |x| in three dimensions, the effective classical potential in Eq. (2.S) 
can be approximated well by variational perturbation theory |l6|,[l8|,[l9| . 
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B. Effective Classical Hamiltonian 



In order to deal wit h H amiltonians like (2.1) which contain a p • A(x)-term, we must generalize the variational 
procedure. Extending (2.8), we define an effective classical Hamiltonian by the phase space path integral 



'/3H ctt (p ,x ) 



z po~ = (2kK) d Iv' D xV D p6(xo - x(r))5(p - p(r)) . 



-A[p,x]/h 



(2.12) 



with the action (2.3) and the measure (| 
phase space integral 



. This allows us to express the partition function as the classically looking 



(2t:K) d 

where po is the temporal average of the momentum: 



exp[-/3ff cff (p ,x )] , 



1 f*f> 

po=p(r)= ^y drp(r) 



(2.13) 



(2.14) 



The fixing of po is done for the same reason as that for Xo, since the classical expectation value (p 2 ) cl is diverging 
linearly with T, just as (x 2 ) cl . 

In the special case of a standard Hamiltonian ( p.6[ ), the e ffective Hamiltonian in Eq. ( (2.13| ) reduces to the effective 
classical potential, since the momentum integral in Eq. (2.12) can then be easily performed, and the resulting restricted 
partition function becomes 



Z po,x = exp _p 



2M J 



(2.15) 



with the local partiti on fu nction Z x ° = exp[— pV c s{^o)] of Eq. (2.8). Thus the complete quantum statistical partition 
function is given by ( [2.13 ) , with an effective classical Hamilton function 



P 2 

ff e ff(P0,X ) = —2- + Veff(xo). 



(2.16) 



As a consequence of the purely quadratic momentum dependence of _ff(p,x) in fl2.6| ), the po-integral in (2.13) can 
be done, thus expressing the quantum statistical partition function as a pur e co nfiguration space integral over the 
Boltzmann factor involving the effective classical potential I4ff(xo), as m Eq. (2.7). 



C. Exact Effective Classical Hamiltonian for an Electron in a Constant Magnetic Field 

The effective classical Hamiltonian for the electron moving in a constant magnetic field can be calculated exactly. 
We consider a magnetic field B = Be z pointing along the positive z-axis. The only nontrivial motion of the electron 
is in the x— y-plane. In symmetric gauge the vector potential is given by 

A(x) = j(-y,x,0). (2-17) 



The choice of the gauge does not affect the partition function since the periodic path integral (2.2) is gauge invariant. 
Ignoring the trivial free particle motion along the z-direction, we may restrict our attention to the two-dimensional 
Hamiltonian 

#(P,x) = %m~ ^c/*(p,x) + \Mul* 2 (2.18) 

with x = (x,y) and p = (p x ,p y ). Here, lo c = eB/Mc is the Landau frequency, and 

k(p,x) = (x x p) 2 = xp y - yp x (2.19) 

the third component of the orbital angular momentum. The partition function of the problem is given by Eq. ( [2.13 ), 
with D — 2. Being interested in an effective classical formulation, we have to calculate the path integral (2.12). First 
we express the <5-function for the averaged momentum as a Fourier integral 
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£(p - p(r)) 
involving an auxiliary source 



( l t \ 
(W eXP ( v -^' P °) eXp 



drv (£) • p(r) 



(2.20) 



(2.21) 



which is constant in time. Substituting the 5- function in Eq. ( [2.12 ) by this source representation, the partition 
function reads 



z po,x = I d 2^ exp ( . pQ ) A V > 2 xV 2 pS ( Xo - x (r)) 



x exp 



dT Hp(t) ' *( T ) + H (P( r )' x ( r )) + v °(£) ' P( T )1 } 



(2.22) 



Evaluating the momentum integrals and utilizing the periodicity property x(0) = x(?i/3), we obtain the configuration 
space path integral 



Z P0,x = Hm / d 2^ exp _ t . p Q 



1 p "'' " (x " ~ x(7,! 



x exp • 



(It 



^x 2 (t) + 1m0 2 x 2 (t) - £mw c (x(t) x x(r)) 2 + x(r) ■ 



(2.23) 



where the source Vo coupled to the momentum in ( 2. 22] ) has turned to a source ji coupled to the path in configuration 
space |2Q], with components 



. / j.\ M iu> c M , . 

Ji(€) = yw c (« 0!( ^),-!) 0l (O) = "2^-(-^ 1 & ! ). 



(2.24) 



We have introduced an additional harmonic oscillator in Eq. fl2.23 ) which will turn out to be useful at intermediate 
stages of the development. At the end of the calculation , onl y the limit £1 — > will be relevant. 
Expressing the 5-function in the path integral of Eq. ( 2.23| ) by the Fourier integral 



(5(x - x(t)) = 



d 2 n 

—- exp [ -k • x ) exp 



1 f h0 

- drhW ■ X ( T ) 



with the new source 



(2.25) 



the partition function ( [2.23| ) can be written as 



z po,x = ^ J d 2^ exp ( . po 



ja(«) = 



M 



d 2 K 



2% 2 p J J (2tt) 



exp -k • x Z n [J (€,«)] 



The functional Zq[J(£, k)] is defined as the configuration space path integral 

Zn[J(€,K)] =y V 2 xexp -- j dr j dr' x(r) G^t. r') x(r') - - j drJ(£,/c)-x 



(r) 



where we have introduced the combined source J(£,k) = ji(£) +J2(k)- Formally, the solution reads 

-a/3 /-ft/3 



Zo[J(€,k)] =^n[0]exp 



-, /-rap ^.n,p 

-—5/ dr / dr'J(|,«)G(T,r')J(^,K) 
2ft Jo Jo 



(2.26) 



(2.27) 



(2.28) 



(2.29) 
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where G(r, r') is the matrix of Green functions obtained by inverting 



G- 1 (r,r') = 



M 



r> 2 - 



' dr 



dr 2 



(2.30) 



The inversion is easily done in frequency space after spectrally decomposing the ^-function into the Matsubara 
frequencies oj m — 2i:m/Ti(3, 



oo 



.(t-t') 



m— — oo 



The result is 



G(u m ) 



M det G V ^ra + fi2 



At this point, the additional oscillator in Eq. (2.23) proves useful: It ensures that the determinant 

dot GK) = (tul + fl 2 ) 2 + w 2 c w 2 m 
is nonzero for m — 0, thus playing the role of an infrared regulator. The Fourier expansion 



1 oo 



u; m )e 



-ICJrn, (T— T ) 



771— — OO 



u T f G xx (t,t') G xy (T,T>) 

'^ T) -{G yx (T,T<) G yv (r,r') 



yields the matrix of Green functions 



which inherits the symmetry properties from the kernel d2.30|) : 

G xx (t, t) = G yy (T, r'), G xy (r, t') = -G vx (t, t'). 

A more detailed description of these Green functions is given in Apps. ^ and [b|. 

Since the current J does not depend on the euclidean time, the expression ( |2.29 ) simplifies therefore to 

r h/3 r hi3 

■J 2 (£.k) I 

h 

The Green function has the Fourier decomposition 



Z n [J(€,K)] - Z n [0]exp 



-jJ^k) / dr \ dr' G xx (t,t') 
n Jo Jo 



-iuj m (r— r ) 



where fi± are the frequencies 



'r> 2 



-W 2 ± W r \ Q 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



(2.35) 



(2.36) 



(2.37) 



(2.38) 



(2.39) 



The ratios in the sum of ( 2.38 ) can be decomposed into two partial fractions, each of them representing a single 
harmonic oscillator with frequency f2+ and f2_, respectively. The analytic form of the periodic Green funct ion o f a 
single harmonic oscillator is well known (see Chap. 3 in and we obtain for the present Green function ( 2.39| ): 



G xx (t,t') = 



1 / h(3 n\ -fl 2 cosh Q + (|r-r'| -hfi/2) h/3 fP - Q 2 coshfi_(|r - r'| - ft/3/2) 



M(3 \2n + n 2 + - n 2 _ 

By writing the determinant ( 2.33| ) as 



sinh h(3VL + /2 



2n_ n 2 



Q 2 



sinh npn_/2 



(2.40) 
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det G(w m ) = (u? rn + n 2 + ){io 2 m + n 2 _) 



(2.41) 



and summing over the logarithms of this, we calculate the partition function as a product of two single harmonic 
oscillators: 



Zr 



Z n [0] 



1 



2 sinh hfin + /2 2 sinh ' 



(2.42) 



The results ( 2.40Q and ( 2.42) ) determine the generating funct ional (2.37). The euclidean time integrations are then 
easily done, and subsequently the k- and ^-integrations in ( 2.27 ). As a result, we obtain the restricted partition 
function 



Z Po.x _ lim ex J a 



1, sinh h(3Vl + /2 sinh hf3tt-/2 p 2 1 , . . 1 , , 1,,„ 2 , 



If we now remove the additional oscillator by taking the limit fi — > 0, we find from ( 2.3S ): f2 + — > oj c , Q_ 
therefore 

sinh7j<9fi + /2 smh?i/3cj c /2 , sinh?i/3fJ_/2 

lllll — — — ; = — ; , lim — — — ; = 1. 



n^o npn+/2 



hPtUa/2 



(2.43) 
0, and 

(2.44) 



Recalling the definition (2.12), we identify the exact effective classical Hamiltonian for an electron in a magnetic field 



1 smhh(3u c /2 pi 1 
gc ff (Po,x ) = -ln h ^ c/2 + — -- 



w c / z (po,x ) + -Mw c x 



(2.45) 



Integrating out the momenta in Eq. (2.13), the configuration space representation ( |2.7| ) for the partition function 
contains the effective classical potential for a charged particle in the plane perpendicular to the direction of a uniform 
magnetic field 



1 sinh h(3uj c /2 
V cS x ) = - In — — — — . 

p nfjuj c /2 

Note that this is a constant potential. 

Denoting the area J cL 2 xq by A, we find the exact quantum statistical partition function 



(2.46) 



Z 



A Hf3uj c /2 
X 2 h sinh?i/3cj c /2' 



(2.47) 



After these preparations, we can turn our attention to the system we want to study in this paper: the hydrogen atom 
in a uniform magnetic field, where the additional Coulomb interaction prevents us from finding an exact solution for 
the effective classical Hamilton function. 



III. HYDROGEN ATOM IN CONSTANT MAGNETIC FIELD 



The zero-temperature properties of the hydrogen atom without external ficls are exactly known. For the quantum 
statistics at finite temperatures, an analytic expression exists, but it is hard to evaluate. It is easier to find an accurate 
approximative result with the help of variational perturbation theory fl9| ] . Similar calculations have been performed 
for the electron-proton pair distribution function which can be interpreted as the unnormalized density matrix Jl8[ . 

Here we extend this method of calculation to the hydrogen atom in a constant magnetic field. This extension 
is quite nontrivial since the weak- and strong-field limits will turn out to exhibit completely different asymptotic 
behaviours. Let us first generalize variational perturbation theory to an electron in a constant magnetic field and 
arbitrary potential. 
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A. Generalized Variational Perturbation Theory 



We consider once more t he effective classical form ( 2.13 ) of the quantum statistical partition function, which requires 
the path integration (2.12) in phase space. Fluctuations parallel and vertical to the magnetic field lines are now both 
nontrivial, and we must deal with the full three-dimensional system and the components of the electron position and 
momentum are now denoted by x = [x, y, z) and p = {p Xl p v , p z ). For the uniform magnetic field pointing along the 
z-axis, the vector potential A(x) is used in the gauge ( 2.17 ). Thus the Hamilton function of an electron in a magnetic 
field and an arbitrary potential V^(x) is 



1 



1 



F(P ' X) = 2M ~ 2^ 2(P ' X) + 8 Mw c 2 x 2 + ^(x). 



(3.1) 



The orbital angular momentum l z (p,x) was introduced in Eq. ( 2.1 9| ) , and the Landau frequency lu c below Eq. (2.18). 
The importance of the separation of the zero frequency components x and po was discussed in Sect. O. Their 
divergence with the temperature T prevents a perturbative treatment. Thus it is essential to set up the perturbation 
theory only fo r th e fluctuations around Xo and po. For this we rewrite the action functional ( |2.3[ ) associated with the 
Hamiltonian (3.1) as 



^l[p,x] =^P- X0 [ P ,x]+Ant[p,x], 



(3.2) 



where we have introduced the fluctuation action 

cft/9 



^• X0 [p,x] = j dr {-»tp(r) - Po] • x(r) + ^[p(t) - Po] 2 + ^h',(pW - Po,x(r) - x ) 
+ Um\ 2 [x» - x^] 2 + \Mn\[z{r) - z ] 2 }, 



in which x^ = {x, y) denotes the transverse part of x. The interaction is now 



/■ft/3 

^ int [p,x]= / drMnt(p(r),x(r)) =^l[p,x]-^' X0 [p,x] 
Jo 

with the interaction potential 

Mnt( P (r),x(r)) = ± {p 2 (r) - [p(r) - Po ] 2 } + \^{r) x x^r) 



(3.3) 



(3.4) 



-io ±1 (p^(r) - p ^) x (x^r) - xjj-) + \m^\t) 
-hml 2 [x» - x^] 2 - \m£1\[z{t) - z ] 2 + V(x(r)), 



(3.5) 



where p 1 - — (p x ,Py)- The frequencies ft = r2j_2, ^||) are for the moment arb itrary. The decomposition (3.2) 

forms the basis for the variational approach, where the first term in the action ( [3.2] ) allows an exact treatment. The 
transverse part was given in Sec. II C and the longitudinal part is trivial, since it is harmonic with frequency fin. The 
associated partition function is given by the path integral 



= d> V^xV'pd^o - x(r))d(p - P(r))e~^ 



[p,x]/ft 



which can be performed. Details are given in Appendix ^|. The result is 



zP0 , X0 = npn + /2 npn_/2 npn l{ /2 

n sinh h/3n+/2 sinh Ti(3Vt- /2 sinh h{3Cl\\ /2 ' 



(3.6) 



(3.7) 



where auxiliary frequencies are composed of the frequencies f2j_i,f2j_2 m t ne action (3.3) as 



n±(o ±1 ,o ±2 ) = -|fi ±1 ±n ±2 | 



(3.8) 



This partition function serves in the subsequent pertubation expansion as trial syst em w hich depends explicitly on 
the frequencies £1. The correlation functions are a straightforward generalization of (2.35) to three dimensions: 
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/G x £(r,r') G x °(t,t') 
G x °(r,r')= G*° x (t,t') G*£(t,t') 

V GS(r,r') 

whose explicit form is derived in App. |^. 

The independent action in Eq. ( |3.2[ ) is treated perturbatively. Writing the partition function ( 2.12| ) as 



(3.9) 



Z^" = (2nh) 3 j>V'\v 3 p5{*o - x(r))5(po - p(r)) exp |-i^P- Xo [p,x]| exp |-i jf*", 



-/r^ nt (p(r),x(T))| , 
(3.10) 



the second exponential is expanded into a Taylor series, yielding 



ZP- X0 = (2ttH) 3 <£v' 3 xV 3 P 6(x -x(T))5(p a -p(T)) exp j-i^^^x] 

I rW 1 /-ft/3 ,hp 

drMnt(p(r),x(r)) + — 2 / dn / dr 2 \4ut(p(n), x(r 1 ))Vint(p(r 2 ), x(r 2 )) 
" Jo 2\n Jo Jo 

Defining harmonic expectation values with respect to the restricted path integral as 

( ■ • • )£°' X ° = |S£ f V'\V 3 V . . . <5(x - xM)<5(Po - W)) exp {-i^°' X ° [P, x] 



the perturbation expansion for the partition function (3.11) reads 

(— l) n / ( r h P 



7 po,xo _ 7Po,x ( 1) /If 



n \ Po,x 



drV- nt (p(r),x(r)) 



(3.11) 



(3.12) 



(3.13) 



This power series expansion can be rewritten in the exponential form 



Z po,x _ z po,x J ( 1) 

J -*n ex P \ 2^ 



n=0 



drMnt(p(r),x(r)) 



(3.14) 



n.c 



where the subscript c on the expectation values indicates cumulants. The lowest cumulants are related to the full 
expectation values as follows: 

(OxkCnMn))^'* - (Oi(p(n),x(n)))r. 

(O 1 ( P (r 1 ),x(r 1 ))O 2 (p(r 2 ),x(r 2 )))P- X0 = ( CMpfa), x(r 1 ))0 2 (p(r 2 ), x(r 2 )) ^ 



;O 1 (p(r 1 ),x(r 1 )))P- X0 (O 2 (p(r 2 ),x(r 2 ))) 



r2 > 



(3.15) 



where Oj(p(rj), x(tj)) deno tes a ny observable depending on position and momentum. Recalling the r elation ( 2.12| ) 
between partition function ( [3.14 ) and effective classical Hamiltonian flcff(Poi x o), we obtain from ( 3.14 ) the effective 
classical Hamiltonian as a cumulant expansion: 



#eff(po,x ) = --j^ z n 



a™n! 



n— 1 



ft /J 



dr Vint (p(r),x(r)) 



n v po.xo 



(3.16) 



Up to now, we did not make any approximation. The expansion on the right-hand side is an exact expression for the 
effective classical Hamiltonian for any fl. 

For systems with a nontrivial interaction, we are capable of calculating only some initial truncated part of the series 
( |3.16 ), say up to the Nth order, leading to the approximate effective classical Hamiltonian 



W[f ) ( Po ,x ) = -iln^- x « 



if HT 

B^> h n n\ 

n—l 



dr VI„t (p(r),x(r)) 



(3.17) 



8 



This depends explicitly on the three parameters ft. Since the exact expression ( [3.16] ) is independent of ft, the 

best approximation for Hf^(po,xo) should depend on fi minimally. Thus the optimal solution will be found by 
determining the parameters from the conditions 

V W#°(po,Xb) = 0. (3.18) 
Let us denote the optimal variational parameters to Nth order by 

nW - (^ ) (Po,x ),^ ) (Po ! x ),r![ | Ar) (po,xo)) . (3.19) 
Inserting these into Eq. ( 3.17j ) yields the optimal effective classical Hamiltonian WW(po,xo). 



B. First-Order Effective Classical Potential 



The first-order approximation of the effective classical Hamiltonian ( 3.17 ) reads 



W&Wxo) = -iln^ Q,X ° - ( V int (p,x))P- X0 . (3.20) 

In writing the last term we have used the fact that, as a consequence of the time translation invariance of the system, 
the first-order expectation value of Vint(x) is independent of the euclidean time r. 

In order to calculate H^(po,xo), we use the two-point correlation functions derived in App. ^, and the vanishing 
of the linear expectations, e.g. 

(p*(t)-Po x Z°'* =0 (3.21) 

to find 

HgWxo) = ^ ~ ~WcWP0,x ) + \M^ c (xl + yl) + Wtf^xo), (3.22) 
where we have collected all terms depending on the variational parameters Q. in the potential 



<(xo) = -iln^°' X0 + (w. - Ou) &i(xo) - j (ni 2 - c c 2 ) ai(x ) - ^faf^o) + (K(x) )P- X ° . (3.23) 



The quantities a^_(xo) and ajj(xo) are the transverse and longitudinal fluctuation widths 

ai(x o ) = GPS' Xo (0), ^(x o )=GP- Xo (0), 6i(x )=GP«; x °(0). (3.24) 



Note that the potential (3.23) is independent of po. This means that the approximation (3.22) to the effective classical 
Hamiltonian contains no coupling of the momentum po to a variational parameter ft, such that the optimal 
determined by minimizing H^\po,'Xo) is independent of po. We may therefore integrate out po in the phase space 
representation of the first-order approximation for the partition function 

.g-^too.xo) (3 25 ) 



„ (1) = f d 3 x d 3 Po p -ff W W (p0|Xo) 
(2irk 3 ) 



to find the pure configuration space integral 



Z V> = f ^°e-^ ! W, (3.26) 



in which Ifi 1 ' (x ) is the first-order approximation to the effective classical potential of an electron in a potential V(x) 
and a uniform magnetic field. 
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C. Application to the Hydrogen Atom in a Magnetic Field 



We now apply the formulas of the preceding section to the Hamiltonian (3.1) with an attracting Coulomb potential 

„2 



V(x) = - 



47T£n X 



(3.27) 



where |x| is the distance between the electro n and the proton. The only nontrivial problem is the calculation of 
the expectation value ( V (x(t)) )q in Eq. (3.23). This is done using the so-called smearing formula, which is a 
Gaussian convolution of V^x). This formula was first derived by Feynman and Kleinert |l7]], and exists now also in an 
extension to arbitrary order JT^JIsJ. The generalization to position and momentum dependent observables was given 
in the phase space formulation pQ]. We briefly rederive the first-order smearing formula. The expectation value is 
defined by 



> '< *< t> ) i : ^r a = SS£ iv'\v 3 pV^{r')) 5(x - x(r))J(p - p(r))e-^ 
i J 



yPo 



(3.28) 



Now we substitute the potential by the expression 



V(x(t'))= / d 3 iy(x)S(x-x(T'))= / d 3 xV(x) 



J d 3 xV(x) J ^-^3 exp[iK- (x-x )]exp|--iy air j(r) • [x(r) - x ] j , 



where we have introduced the source 

j(r) = ifiK8{T — t') 
Inserting the expression ( ^29| ) into Eq. ( ^28|) we obtain 

(27T) 

with the harmonic generating functional 



(y(x(r'))>P°- X ° 



7 p 
z f2 



^ / d 3 xV(x) J exp^.(x-x )] ^ x °[j], 



ZP- X0 [j] = (2nh) 3 fv' 3 xV 3 P 5(* - x(r))5( P0 - p(r)) exp J -i^' xo [p,x] - \ dri(r) ■ [x(r) - x ] 



The solution is 



^°' Xo [j] = ^°' Xo exp 



-S dr dr'j(r)G x o(r,r')j(r') 
Ztl Jo Jo 



(3.29) 
(3.30) 

(3.31) 
(3.32) 
(3.33) 



with the 3 x 3- matrix of Green functions of Eq. (3.E). The properties of th e Green functions are discussed in the 
Appendices |X] and |b| Expressing the source j(r) in terms of k via Eq. ( 3.30 ) and performing the r-integrations, we 
arrive at 



V(r(r'))) P n °' X0 = J d 3 xV(x) J exp{iK. [x-x ]} exp 



-kG Xo (0)k; 



(3.34) 



Recognizing that G x £(0) = G x °(0) vanish, the /«-integral is easily calculated and leads to the first-order smearing 
formula for an arbitrary position dependent potential 



l\\ \ Po,x 



:vw))>n 



(2tt)3/2 a ^(x ),/a2(x 



d 3 x V(x) exp 



(x - x ) 2 + (y ~ yo) 2 (z - z Q ) 2 



2ai(x ) 



2af(xo) 



(3.35) 



the right-hand side containing the Gaussian fluctuation widths (3.24). 
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For the Coulomb potential (3.27) that we are interested in, the integral in the smearing formula ( |3 .35 ) can not be 



done exactly. An integral representation for a simple numerical treatment is 



Po,x 



47r£ o |x| / n 47re V 7r II J a|(x ) + $ 2 [a^(x ) - a|(x )] 




X { L|(x )+e 2 K(xo)-a|(xo)] + ^(xoT J J ■ (3 - 36) 



With this expression we know the entire first-order effective classical potential ( 3.23| ) for an electron in a Coulomb 



potential and a uniform magnetic field which has to be optimized in the variational parameters ft. 

IV. RESULTS 

We are now going to optimize the effective classical potential by extremizing it in ft at different temperatures and 
magnetic field strengths. In the zero-temperature limit this will produce the ground state energy. 

A. Effective Classical Potential for Different Temperatures and Magnetic Field Strengths 

The optimization of (xq) proceeds by minimization in ft and must be done for each value of xq. Reinserting 



the optimal parameters nW(xo) into the expressions ( 3.25 ) and ( 3.36| ), we obtain the optimal first-order effective 



classical potential W^'feo). The calculations are done numerically, where we used natural units h = e 2 /Aite^ = 
c = M = 1. This means that energies are measured in units of eo = Me 4 / (47T£o) 2 ?i 2 = 2 Ryd w 27.21 eV, temperatures 
in eo/ks ~ 3.16 x 10 5 K, distances in Bohr radii as — (4:Trso) 2 h 2 /Me 2 ~ 0.53 x 10~ 10 m, and magnetic field strengths 
in Bq = e 3 M 2 /h 3 (4:TTeo) 2 w 2.35 x 10 5 T = 2.35 x 10 9 G. Figure |l| shows the resulting curves for various magnetic 
field strengths B and an inverse tempature j3 = 1/T = 1. Examples of the lower temperature behaviour are shown 
in Fig. H for = 100. To see the expected anisotropy of the curves in the magnetic field direction and in the plane 
perpendicular to it, we plot simultanously the curves for W^(no) transversal to the magnetic field as a function of 
po = v^o + Va a t z = (solid curves) and parallel as a function of z Q at p — (dashed curves). The curves become 
strongly anisotropic for low temperatures and increasing field strengths (Fig. ||). At a given field strength B, the two 
curves converge for large distances from the origin, where the proton resides, to the same constant depending on B. 
This is due to the decreasing influence of the Coulomb interaction which shows the classical 1/r-behaviour in each 
direction. When approaching the classical high-temperature limit, the effect of anisotropy becomes less important 
since the violent thermal fluctuations do not hav e a preferred direction (see Fig. pi). For pa — > oo or zq — > oo, the 



expectation value of the Coulomb potential ( 3.36 ) tends to zero. The remaining effective classical potential 



W^xo) — > -ilnZP- X0 + ( Wc - n xl ) b\ - I (Q 2 ±2 - lo 2 ) a\ - -Mfifoj} (4.1) 

is a constant with regard to the position Xo, and the optimization yields fl^l = = u c and fl^ = 0, leading to 
the asymptotic constant value 

The B = -curves are of course identical with those obtained from variational perturbation theory for the hydrogen 
atom ATI. 
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FIG. 1. Effective classicai potentiaf as a function of the coordinate po = y/%o + v\ perpendicular to the field lines at Zq = 
(solid curves), and parallel to the magnetic field as a function of Zo at po = (dashed curves). The inverse temperature is 
fixed at /3 = 1, and the strengths of the magnetic field B are varied (all in natural units). The small figure enlarges the range 
< po , zo < 1 with noticeable anisotropy. 
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FIG. 2. Analogous plot to Fig. [j], but at the larger inverse temperature (5 = 100 
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B. Ground State Energy of the Hydrogen Atom in Uniform Magnetic Field 

In what follows we investigate the zero-temperature behaviour of the theory. Figures |l| and || show that the minimum 
of each potential curve lies at the origin. This means that the first-order approximation to the ground state energy 
for a fixed magnitude of the magnetic field B is found by considering the zero-temperature limit of the first-order 
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effective classical potential in the origin 



= lim W {1) (0). 



Thus we obtain from Eq. ( 3.23] ) the variational expression for the ground state energy: 

o 



4ft 



ft 



an ii 



Atteq \ x 



1 



(4.3) 



(4.4) 



where the expectation value for the Coulomb potential (3.36) can now be calculated exactly since the exponential in 
the integral simplifies to unity: 



ft||ft±2 / 20 1| 

■ arctan - 



fti — ft 



ft 



On, 



ft||ft±2 , 1 + 2 v /2ri||/S2j_ 2 — 1 

In ■ 



2i y fin — n_L2 i -i^/2n ll /n_ L2 - 1 



2ft|, > n_L2, 

2ft|| = ft^2, 
2ft|| < ftj_2- 



(4.5) 



The equations (4.4) and (4.5) are independent of the frequency parameter ftj_i such that the optimization of the 
first-order expression for the ground state energy (4.4) requires the satisfying of the equations 



aft 



±2 



90,1 



(4.6) 



Reinserting the resulting values ft^ an d into Eq. (4.4) yields the first-order approximation for the ground state 



(i) 

±2 



L0 C 



energy E^(B). In the absence of the Coulomb interaction the optimization with respect to ft^2 yields ft 
rendering the ground state energy E^ X \B) = u) c /2, which is the zeroth Landau level. An optimal value for ft|| does 
not exist since the dependece of the ground state energy of this parameter is linear in Eq. (4.4) in this special case. To 
obtain the lowest energy, this parameter can be set to zero (all optimal frequency parame ters used in the optimization 
procedure turn out to be nonnegative) . For a vanishing magnetic field, B = 0, Eq. ( |l.4| ) exactly reproduces the 
first-order variational result for the ground state energy of the hydrogen atom, E^\B = 0) « -0.42 [2Ryd], obtained 
in Ref. @. 

To investigate the asymptotics in the strong-field limit B — > oo, it is useful to extract the leading term uj c /2. Thus 
we define the binding energy 



E(B) 



(4.7) 



which possesses an characteristic strong-field behaviour to be discussed in detail subsequently. The result is shown 
in Fig. H as a function of the magnitude of the magnetic field B, where it is compared with the high-accuracy 
results of Ref. Q. As a first-order approximation, this result is satisfactory. It is of the same quality like other 
first-order results, for example those from the operator optimization method in first order of Ref. B. The advantage 
of variational perturbation theory is that it yields good results over the complete r ang e of the coupling strength, 
here the magnetic field. Moreover, as a consequence of the exponential convergence [|16|, Chap. 5], higher orders of 
variational perturbation theory push the approximative result of any quantity very rapidly towards the exact value. 
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1000 2000 
B/B 

FIG. 3. First-order variational result for the binding energy as a function of the strength of the magnetic field. The dots 
indicate the values of Ref. JlJ. The dashed curve shows the simple estimate of Landau-Lifschitz || 0.51n 2 _B, which is closely 
related to the ground state energy of the one- dimensional hydrogen atom 



1. The Weak-Field Case 



We investigate now the weak-field behaviour of our theory starting from the expression (4.7) and the expectation 
value of the Coulomb potential (4.5) in natural units: 



*W(B) = ---[! + -) - — —h(rj) 

»7> m ' 9 A \ 9 J AO V 9tt w; 



B (I 



B 2 



lr)(l. 



40 



2vr 



with 



1 



:hl 



In comparison with Eq. (4.4) we introduced new variational parameters 

2flu 



(I 



n = n ±2 



(4.8) 



(4.9) 



(4.10) 



and utilized, as the calculations for the binding energy showed, that always T) < 1. Performing the derivatives with 
respect to these variational parameters and setting them zero yields conditional equations which can be written after 
some manipulations as 



O 



(i 



1 



y 2-kt] 1 — rj 
1 7] B 2 1 



1 



1 1 



1 - VT^Tj 

1 



0. 



i] 



4 8 Afl 2 ' 2 V 2tt^ y/T^'" 1 + ^/T=rj 
Expanding the variational parameters into perturbation series of the square magnetic field B 2 , 



(4.11) 



V 



(B) = J2 VnB 271 , (1(B) = VnB 2 



(4.12) 



n=0 



n=0 
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and inserting these expansions into the self-consistency conditions (4.11) and (4.11) we obta in o rder by order the 
coefficients given in Table |[ Inserting these values into the expression for the binding energy ( |4.8| ) and expand with 
respect to B 2 , we obtain the perturbation series 



-£ e >vB 2n - (4-13) 

n=0 

The first coefficients are also given in Table |. We find thus the important result that the first-order variational 
perturbation solution possesses a perturbative behaviour with respect to the square magnetic field strength B 2 in the 
weak-field limit thus yielding the correct asymptotics. The coefficients differ in higher order from the exact ones but 
are improved in higher orders of the variational perturbation theory ||l6| , Chap. 5]. 



TABLE I. Perturbation coefficients up to order B 6 for the weak-field expansions of the variational parameters and the 
binding energy in comparison to the exact ones of Ref . [[| . 



n 









1 


2 




3 








1.0 


405vr 2 
7168 


- w -0.5576 


1682 8 9657T 4 
1 25881 5488 


w 1.3023 


38869993320757r b 
884272562962432 


w -4.2260 




32 
9tt 


» 1.1318 


99vr 
224 


ss 1.3885 


1293975tt j 
1 9668992 


a -2.03982 


524431667187tt s 
27633517592576 


« 5.8077 


En 


4 

~3^ 


» -0.4244 


9tt 
128 


ss 0.2209 


80197^ 


i -0.1355 


2564498077r b 
322256764928 


0.2435 


e n | 






-0.5 




0.25 


1 835j)08 
_ 192 ~ " 


-0.2760 


« 1.2112 

4608 
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2. Asymptotical Behaviour in the Strong-Field Regime 



In the discussion of the pure magnetic field below Eq. (|^) we have mentioned that the variational calculation for 
the ground state energy which is thus associated with the zeroth Landau level yields a frequency Oj_2 oc B while 



o 



0. Therefore we use the assumption (with Oj_ = 0^2) 

0_l>20||, fi||<B 



(4.14) 



for the consideration of the ground state energy (4.4) of the hydrogen atom in a strong magnetic fi eld. In a first 
step we expand the last expression of the expectation value ( |4.5| ) which corresponds t o th e condition ( 4.14 ) in terms 
of 2f2||/f2_L and reinsert this expansion in the equation of the ground state energy (4.4). Then we omit all terms 
proportional to C/Oj_ where C stands for any expression with a value much smaller than the field strength B. In 
natural units, we thus obtain the strong-field approximation for the first-order binding energy ( [4.7] ) 



2 4 40 1 4 V 7T 20 



(4.15) 



As usual, we consider the zeros of the derivatives with respect to the variational parameters 



de 



(i) 
dO I, 



de, 



(i) 



ao 



which lead to the self-consistence equations 



0|| = — 7=(lnO|| -lnO_L + 2-ln2) 



± = 2\ + B\ 1 



0|| 
[ ^B~2 



(4.16) 

(4.17) 
(4.18) 



Let us first consider the last equation. Utilizing the second of the conditions (4.14) we expand the second root around 
unity yielding the expression 



Oi =B 











2 !L _ 4 ^_ 

ttB tt 2 B 3 



(4.19) 



where the terms are sorted with regard to their contribution starting with the biggest. Since we are interested in the 
strong B limit, we can obviously neglect terms suppressed by powers of l/B. Thus we only consider the following 
terms for the moment: 







/On 
V 7T 



(4.20) 



Inserting this into the other condition (4.17), expanding the corresponding logarithm, and, once more, neglecting 
terms of order 1/-B, we find 



InB-InO 



(i) 



In 2 



(4.21) 



To obtain a tractable approximation for Oil, we perform some iterations starting from 

2 







(i) 



=ln2Be- 2 



(4.22) 



Reinserting this on the right-hand side of Eq. (4.21), one obtains the second iteration ^/Oy 2 ' ) . We stop this procedure 
after an additional reinsertion which yields 







(3) 



ln2Se" 2 - 21n 



ln2Be _2 - 21n 



zln2.Be- 2 



(4.23) 
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The reader may convince himself that this iteration procedure indeed converg es. F or a subsequent systematical 
extraction of terms essentially contributing to the binding energy, the expression ( 4.23 ) is not satisfactory. Therefore 
it is better to separate the leading term in the curly brackets and expand the logarithm of the remainder. Then this 
proceeding is applied to the expression in the angular brackets and so on. Neglecting terms of order ln~ 3 B, we obtain 



d 3) m ^=(\n2Be- 2 



In 21nln 2Be 

4 



(4.24) 



The double-logarithmic term can be expanded in a similar way as described above: 



lnln 2BeT 2 = In 



lnB 1 



In 2 



lnB 



= lnln B 



Thus the expression (4.24) may be rewritten as 

2 



with abbreviations 



n 



(3) 
I 



In B - 21nln B 



In 2 - 2 1 (In 2 - 2) 



lnB 



2a | a 2 
h^B + hFB 



In B 



0(kT 3 B). 



b ] + e>(hi~ 3 B) 



(4.25) 



(4.26) 



a = 2- In 2 sa 1.307, 



6 = ln--2 



-1.548. 



(4.27) 



The first observation is that the variational parameter On is always much smaller than f2j_ in the high B-field limit. 
Thus we can further simplify the approximation (4.20) by replacing 




(4.28) 



without affecting the following expression for the binding energy. Inserting the solutions ( 4.26| ) and (4.28) into the 
equation for the binding energy (4.15) and expanding the logarithmic term once more as described, we find up to the 



order In B: 



■«(B) 



ln 2 B - 4 In B lnln B + 4 ln 2 ln B - 46 lnln B + 2(6 + 2) In B + b 2 



-e>(kr 2 B) 



i 



lnB 



[8 ln 2 lnB 



>lnlnB + 26 2 ] 

(4.29) 



Note that the prefactor 1/ir of the leading ln 2 B-term differs from a value 1/2 obtained by Landau and Lifschitz [||. 
Our different value is a consequence of using a harmonic trial system. The calculation of higher orders in variational 
perturbation theory would improve the value of the prefactor. 

At a magnetic field strength B = 10 5 Bo, which corresponds to 2.35 x 10 10 T = 2.35 x 10 14 G, the contribution from 
the first six terms is 22.87 [2 Ryd]. The next three terms suppressed by a factor ln _1 B contribute —2.29 [2 Ryd], while 
an estimate for the ln _2 B-terms yields nearly —0.3 [2 Ryd]. Thus we find 



e (1) (10 5 ) = 20.58 ±0.3 [2 Ryd]. 



(4.30) 



This is in very good agreement with the value 20.60 [2 Ryd] obtained from the full treatment described in Sec. [V B| . 

Table [H] lists the values of the first six terms of Eq. ( 4.29| ). This shows in particular the significance of the second- 
leading term — (4/7r)lnB mini?, which is of the same order of the leading term (l/7r)ln 2 B but with an opposite sign. 
In Fig. ||, we have plotted the expression 



e L (B) = -\n 2 B 



(4.31) 



from Landau and Lifschitz Q to illustrate that it gives far too large binding energies even at very large magnetic 
fields, e.g. at 2000B cx 10 12 G. 

This strength of magnetic field appears on surfaces of neutron stars (10 10 — 10 12 G). A recently discovered new 
type of neuton star is the so-called magnetar. In these, charged particles such as protons and electrons produced by 
decaying neutrons give rise to the giant magnetic field of 10 15 G. Magnetic fields of white dwarfs reach only up to 
10 6 — 10 s G. All these magnetic field strengths are far from realization in experiments. The strongest magnetic fields 
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ever produced in a laboratory were only of the order 10 5 G, an order of magnitude larger than the fields in sun spots 
which reach about 0.4 x 10 4 G. Recall, for comparis on, th at the earth's magnetic field has the small value of 0.6 G. 

As we see in Fig. |[ the nonleading terms in Eq. (4.29) give important contributions to the asymptotic behaviour 
even at such large magnetic fields. It is an un usual property of the asympto tic be haviour that the absolute value of 
the difference between the Landau-expression (4.31) and our approximation (4.29) diverges with increasing magnetic 
field strengths B, only the relative difference decreases. 



TABLE II. Example for the competing leading six terms in Eq. (kL29) at B = 10 5 B » 2.35 x 10 14 G 



(l/7r)ln 2 B 


-(4/7r)lnB lnlnB 


(4/tt) ln 2 lnB 


-(4&/7r)lnln£ 


[2(6 + 2) /tt] InB 


b 2 /n 


42.1912 


-35.8181 


7.6019 


4.8173 


3.3098 


0.7632 
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V. SUMMARY 



We have calculated the effective classical potential for the hydrogen atom in a magnetic field. For this we have 
generalized variational perturbation theory to make it applicable to physical systems with uniform external magnetic 
field. 

The effective classical potential containing the complete quantum statistical information of the system was deter- 
mined in first-order variational perturbation theory. For zero-temperature, it gave the energy of the system. Our 
result consists of a single analytic expression which is quite accurate at all temperatures and magnetic field strengths. 
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APPENDIX A: GENERATING FUNCTIONAL FOR PARTICLE IN MAGNETIC FIELD AND 

HARMONIC OSCILLATOR POTENTIAL 

For the determination of the correlation functions of a system, we need to know the solution of the two-dimensional 
generating functional in the presence of an external source j = {j x ,jy)'- 



-A*° [xtJ]/R 



The action of a particle in a magnetic field in z-direction and a harmonic oscillator reads 



^ X0 [x;j] 



dr 



^x 2 (r) - ijfu([x(T) - xo] x x(r)) z + ^M0 2 [x(r) - x ] 2 + j(r) • (x(r) - x ) 



(Al) 



(A2) 



The position dependent terms are centered around x = (%q, yo), which is the temporal average of the path x(r), and 
thus equal to the zero frequency component of the Fourier path, is 



x(r) = xo + ( x ^ mT + <e-^ T ) 



(A3) 



with the Matsubara frequencies uj m — 2-Km/fif3 and complex Fourier coefficients x m = xjjj + ix™. Introducing a 
similar Fourier decomposition for the current j (r) with Fourier components j m and using the orthonormality relation 



1 C hf3 
— / dre 1 ^-""' 7 = S„ 



W Jo 

the generating functional can be written as 

dx™dxZdy™dy. 



z x °\i]= n 

ra— 1 



nn 

Hi > ^— A-m (x m ^ m )/H 



(A4) 



(A5) 



with 



+2hf3(x m j Xm + x m j Xm + y m jy m + Dm 3y m )■ 



(A6) 



Expression (A5) is equivalent to the path integral (Al) and we obtain after performing the integrations and retrans- 
forming the currents 
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1 fW 
* m = WJo dTiiT)e 



—ZUJm T 



the resulting generating functional 



Z x «[j] = Z*° exp S ^2 y o dr J dr> j(r)G*° (t,t>)j(t') \ 



with the partition function 

oo 

z x ° = z x °[o] = rj 

m— ] 

and the 2 x 2-matrix of Green functions 



+ + ^2)2 



GS(r,r') G x °(t,t') 
G x S(r,r') G x £(t,t') 



(A7) 



(A8) 



(A9) 



(A10) 



The elements of this matrix are position-position correlation functions what can be easily proved by applying two 
functional derivatives with respect to the desired component of the current to the functional (Al), for example 



GZ° x (t,t') = ((x(t)-x )(x(t')-x ))> 



where we have defined expectation values by 



\2 
A, 



Z^[S\5j x (T)5j x {T>) 



j=o 



( . . . ) X0 = ^ j>&* . . . *(xo - x(r)) e -W^ 



(All) 



(A12) 



From the above calculation we find the following expressions for the Green functions in Fourier space (0 < r, r' < hp) : 

G x £(t,t') = <£(r)i(r')) x ° = G x «(r,r') = ( y(r) y(r') ) x ° 



E 



-lOJ-m (t — T ) 



G x »(t,t') = (i{r)y{r')r = -G?°(t,t') = - ( y(r) x(r') ) xo 



2cj 



E 



-iUJm { T ~ r ' ) 



(A13) 



(A14) 



where, for simplicity, x(r) = x(r) — xo. It is desirable to find analytical expressions for the Green functions and the 
partition function (A9). All these quantities possess the same dominator which can be decomposed as 



co 2 Lo 2 m + (cj 2 m + n 2 ) 2 = (co 2 m + n 2 + )(u 2 n + n 2 _) 



with frequencies 



+ ^Lo 2 ±uj^n 2 + 



(A15) 



(A16) 



Therefore the partition function (A9) can be split into two products, each of which known from the harmonic oscil- 
lator U, Chap. 5]: 



Z x ° = 



n 



coi + n 2 + 



n 



u 2 +W 



npn + /2 npvi-/2 



smhhpn + /2 smhnpn^/2' 



(A17) 



Now we apply the property ( A15 ) to decompose the Green functions ( |A13| ) into partial fractions, yielding 
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1 / 00 1 
G^( T , t') = G^S(t, t') = -— ai V -= -r e -' l ^ T - T ) 

\ m = — OO "* ~r 



«2 



E 



j-iwmfr-r') 



m— — oo 



o 2 



with coefficients 



Ql 



n\-n 2 
n 2 -n 2 _ 



(A18) 



(A19) 



Following Ref. |T1| Chap. 3], sums of the kind occuring in expression ( |A16| ) are spectral decompositions of the 
correlation function for the harmonic oscillator and can be summed up: 



E 



i 



4 + & ± 



h[3 coshfi ± (|r-r'| -h/3/2) 



2fl± smhh(3n±/2 
Thus, the xx- and yy-correlation functions can be expressed by 

i / h/3 fi% -n 2 cosh^+(|r-T'| -a/3/2) h/3 n 2 _-n 2 coshn-(\T-T'\- h/3/2) 1 



(A20) 



M(3 \2Vt + 02 



n 2 



sinh h/3n + /2 



20_ ni - n 2 _ 



sinh npn^/2 



n 2 



(A21) 



where, from Eq. (A16), Q± — il±(u!,Q) are functions of the original freque ncies u> from the magnetic field and O 
from the additional harmonic oscillator ([A2]). It is obvious that expression (A.21) reduces to the Green function of 
the harmonic oscillator for oj — ► 0: 



^^=M^\ 2 



coshOflr — r' 



K/3/2) 



sinh h/3n/2 



- 1 



(A22) 



with i £ {x,y\. In this limit, the partition function (A.17) turns out to be the usual one |T^, Chap. 5] for such a 
harmonic oscillator 



u^q sinh hpn/2' 



(A23) 



It is worth mentioning that with the last term in Green function ( A21 ) the classical harmonic fluctuation width 

cl 1 



GC 



2\ cl 



(x 2 ) 



M/3fl 2 



(A24) 



is subtracted. This is the con sequence of the exclusion of the zero frequency mode of the Fourier path (|Aj) in 
the gen eratin g functional ( ]Al| ). The necessity to do this has already been discussed in Sect. |l[ The other terms 
in Eq. (A21) are those which we would have obtained without separation of the cco-component. Thus these terms 
represent the quantum mechanical Green function containing all quantum as well as thermal fluctuations. It is a nice 
property of all Green functions discussed in this paper that 



(A25) 



Such a relation exists for all other Green functions appropriately, including momentum-position correlations which 
we consider subsequently. 

The knowledge of relation ( A2C| ) make s it quite easy to determine the algebraic expression for the mixed xy- 
correlation functions. Rewriting Eq. ( A. 14) as 



G^t>) = -G^r') = — 



1 







M/3 Q 2 , - Q 2 _ dr 



E 



Q 2 + 



E 



»(t-t') 



f! 2 



(A26) 



and applying the derivative with respect to t to relation ( A2C| ), we obtain the following expression for the mixed 
Green function: 



Gxy{ T ' T ') — ~ CyS( r : T ') 

hw l 



2iM n 2 + - n 2 



{6(r - r')[g n+ (r, r') - g n _ (r, r')] - 6(r' - r)[g n+ (r', r) - g Q _ (r', r)} } 



(A27) 
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where we have used the abbreviation 



, sinhfi±(r-r'- %0/2) 
9n± ^~ sinh?i/?r! ± /2 ' 

Note that classically (xy) cl = such that Eq. ( |A25| ) reduces to 

G? y (T,T') = G^(T : T<). 



t,t' e (o,n/3). 



The Heaviside function in Eq. ( A27 ) is defined symmetrically: 

e(r-r') = 



1 t>t', 
1/2 r = r', 







t < r 



In the quantum mechanical limit of zero-temperature (/3 — > oo), the Green function (A21) simplifies to 



/ 1 SI 2 SI 2 

lim G^(r,r') = lim G*°Jt,t') = — — ^ 9 



^ — Oo- It— 7 



1 £1 2 - f7 2 

o! n 2 + -n 2 _ 



,-n_|r-7 



while in Eq. ( |A27| ) only gn ± (r, r') changes: 



lim g n± (t, t ) = -e 



'1 = _ p -"±(T-r') 



(A28) 



(A29) 



(A30) 



(A31) 



(A32) 



APPENDIX B: PROPERTIES OF GREEN FUNCTIONS 



In this section we list properties of the Green functions (A21) and (A27) which are important for the forthcoming 
consideration of the generating functional with sources coupling linearily to position or momentum in Appendix KJ 
For all relations we suppose that < r, t' < %j3. 



1. General Properties 

A first observation is the temporal translational invariance of the Green functions: 

where each of the indices i,j stands for x or y, respectively. For equal times we find 



vi \ 1 / nB n 2 - n 2 , m n 2 - n 2 

Gil(r, T ) = TTia I o*T~ o2 + 02 coth^n+72 



M/3 \2fl+ n 2 - n 2 



20_ n 2 -n 2 _ C0thh ^ 2 -^ 



Moreover we read off the following symmetries from the expressions ( A21 ) and ( A27 ) 



Otherwise, 



G~°(t,tI) = G*°(t',t). 



1 i = 3, 
i^j. 



(Bl) 



(B2) 



(B3) 



(B4) 



Throughout the paper we always use periodic paths. Hence it is obvious that all Green functions are periodic, too: 

G*P(0,t') = G%(hf3,T'), G*°(r,0) = O${T,h0). (B5) 
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2. Derivatives of Green Functions 



We now proceed with derivatives of the Green functions (A21) and (A.27), since these are essential for the derivation 
of the generating functional of position and momentum dependent correlations in the forthcoming Appendix 

Before considering the concrete expressions we introduce a new symbol indicating uniquely to which argument the 
derivative is applied. A dot on the left-hand side means to perform the derivative with respect to the first argument 
and the dot on the right-hand side indicates that to differentiate with respect to the other argument. Having a dot 
on both sides the Green function is derived with respect to both arguments: 



8t 



d 2 G* (r,r') 
drdr' 



Applying such derivatives to the Green functions (A21), we obtain (i £ {x,y}): 

h 1 



•G*°(t,t') = 



2M n\ - n 2 



[6(r - r')g(r, t<) - Q(r' - r)g(r\ r)] = -G'*°(t, t') 



with 



(B6) 



(B7) 



9 (t,t') = (Q 2 + - n 2 )g n+ (T,T') - (Q 2 _ - n 2 )g (l _ (r, r'), 
where g± (r, t) was defined in Eq. ( |A28|) . Performing the derivatives to both arguments leads to the expression 



(B8) 



•G^(t 1 t>) = -G-*:(t 1 t>) + -8(t-tI) 



where we have introduced the partial function 



•a?(ry>~£ 



n 2 + -fl 2 smhO + (|r-r'| - h/3/2) Q 2 _ - fl 2 sinhfi_(|T - r'| - h/3/2) 

^+FT2 — ^-r-^z — " 



sinh hpn+/2 



n 2 



sinh h(3VL-/2 



(B9) 



(BIO) 



which is finite for equal times. 

Applying derivatives with respect to the first respective second argument to the mixed correlation function ( A.27] ) , 
we find: 



°G*° v (t, t) 
and 



Hu) 



1 



2iM ni 



n 2 



coshrj+(|T-r'|- h/3/2) coshfi_(|r-T'| -h/3/2) 



sinh hptt + /2 



sinh hpn_/2 



'G v x {t,t) — —*G X y(r,T). 
Differentiating each argument of the mixed Green function results in 



-G'2(t,t>) (Bll) 



(B12) 



'G'*°(t,t') = 



ihuj 



2M n 2 + - n 2 _ 



[9(r - r')h(r, r') - 6(r' - r)h(r',r)} = -'G'* x (t, t') 



with 



h(r, t') = n 2 +9n+ (r, t') - n 2 _g n _ (r, r'). 



An additional property we read off from Eqs. (B7) and (Bll) is 6 {x,y\)\ 

-1 



•G*°(t,tI) = 'G*<>(t\t)x 
G-r/(r,r') = G-^(r',r)x 



1 



1 



i = 3, 
i + 3, 

i =3; 

i +3- 



(B13) 

(B14) 

(B15) 
(B16) 



The double-sided derivatives flB9|), ( |B10|) , and ( |B13|) imply 



23 



i =3, 



(B17) 



The derivatives (p7|) , ( pTp| ) , ( ^TT] ), and (^T|) are periodic 



•Gf/(r,0) = 

G**°M) = 

•G"«(r,0) 
•G«P(r,0) 



Gt/(r, hp), 
= 'G-*°(r,hP), 

= 'G'*°(r,hp), 



•c^(oy) 

GT/(0,r') 
•G^ o (0,r') 



(0,r') 



(B18) 
(B19) 

(B20) 
(B21) 



APPENDIX C: GENERATING FUNCTIONAL FOR POSITION- AND MOMENTUM-DEPENDENT 

CORRELATION FUNCTIONS 



With the discussion of the generating functional for position-dependent correlation functions and, in particular, 
the Green functions in Appendix ^ and their properties in Appendix we have layed the foundation to derive the 
generating functional for correlation functions depending on both, position and momentum. Following the framework 
presented in an earlier work ]2C| ] , such a functional involving sources coupled to the momentum can always be reduced 
to one containing position-coupled sources only. 

We start from the three-dimensional effective classical representation for the generating functional 



d 3 x d 3 p vPo , Xor . , 
(2^)3 ^ [J ' VJ 



(CI) 



with zero frequency components x = (^oi ^0) = const, and p = {PxQ,Py ,Py ) = const, of the Fourier path 
separated. The reduced functional is 



Z*°' xo \j,v] = (2^) 3 ^P' 3 ^ 3 p ( 5(xo-x(r)) ( 5(po-p(r)) exp j-±4£ 0,X0 [p,x;j,- 



(C2) 



where the path integral measure is that define d in Eq. (|2.4|). Extending the action (2.3) by source terms, considering 
a more general Hamilton function than ( [2.17| ), and introducing an additional harmonic oscillator in z-direction, the 
action functional in Eq. (C2) shall read 



r h/3 rill I 
^' X0 [p,x;j,v] = J dr {-ip(r) ■ x(r) + ^P 2 (r) ^ ±1 l z (p,i) + ~M& ±2 [Z\r) + f(r)] + -MQp 2 (r) 

+j(r)-i(r)+v(r)-p(r)} (C3) 

with shifted positions and momenta 

x = x(r)-x , P = p(r)-p . (C4) 



The orbital angular momentum / z (p,x) is defined in Eq. ( [2.19| ) and is used in Eq. (C3) with the shifted phase space 
coordinates (C4). We have introduced three different frequencies in (C3), £1 = (fi_u, ^X2> where the first both 
components are used in regard to the oscillations in the plane perpendicular to the direction of the magnetic field 
which shall be considered here to point into z-direction. The last component, Qy, is the frequency of a trial oscillator 
parallel to the field lines. 

Due to the periodicity of the paths, we suppose that the sources might also be periodic: 

j(0)=j(ft/3), v(Q) = v(7i/?). (C5) 

Since we want to simplify expression ( |C2| ) such that we can use the results obtained in Appendix [A], the momentum 
path integral is solved in the following. In a first step we reexpress the momentum (5-function in (|C2[) by 



S(P0 - p(r)) = J ^ exp j-~ ^ dr v . [p(r) - p ] 



(C6) 
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where 



(C7) 



is an additional current which is coupled to the momentum and is constant in time. Defining the sum of all sources 
coupled to the momentum by 



V(€,r)=v(r)+v (C), 



(C8) 



the functional ( |C2| ) can be written as 



7P0,X [ 



[],v]= I d^<£V' 3 xV 3 p5(xo-x(T))cxp{ - j I dr - ip(r) ■ x(r) + - ~0^( P (t), x(r)) 



lit) 



P 2 (r) 1 



2M 



+ \mQ 2 ±2 {£ 2 (t) + y 2 (r)} + \mQ\ ~z 2 (t) + j(r) ■ x(r) + V(£, r) ■ p(r) 



(C9) 



where we have used the translation invariance p — ► p of the path integral. To solve the momentum path integral, it is 
useful to express it in its discretized form. Performing quadratic completions such that the momentum path integral 
separates into an infinite product of simple Gaussian integrals which are easily calculated, the remaining functional 
is reduced to the configuration space path integral 



y^exp M.J drV^for) y^<5(x -x(r))exp|-^P- Xo [x;j,V]| 



(CIO) 



with the measure ( 2.10] ) for D = 3. The action functional is 



^P- X0 [x;j,V]- / dr 
Jo 



^x 2 (r) + iiMfixi {i(r)y{r) - y{r)i{r)} + l -M (tf ±2 - tf ±1 ) {x\r) + f{r)} 



+-Mnp 2 { T ) + x (T) 



j x (T) + ±Mn ±1 v y (s,T) 



+y(r) 



j v (r)--MSlxMi,T) 



iM r n/3 

l — J drx(r).V(€,r), 



where the last term simplifies by the following consideration. A partial integration of this term yields 



fit 



drx(r) • V(C, r) = - / dr (x(r) - x ) • V(£, r). 



(Cll) 



(C12) 



The surface term vanishes as a consequence of the periodicity of the path and th e sou rce. This periodicity is also the 
reason why we could shift x(r) by the constant x on the right-hand side of Eq. ( C12| ). Obviously, the importance of 
this expression lies in the coupling of t he ti me derivative of V(£, r) to the path x(t). Thus, V(£, t) can be handled 
like a j(r)-current an d the action ( Cll ) can be written as 



A%>*° [x; j, V] = A%>^> [x; J, 0] = Ag'* [x; 0, 0] - \ dr x(r) • J(£, r) 



(C13) 



with the new current vector J(£,t) which has the components 



■/«(€, t) = Mr) + 2 M n±iV y (S, t) - iMV x {Z, t), 
M& r) = Mr) ~ \Mn ±1 V x (t t) - iMV y {£, r), 



(C14) 



Mi,r) = j z {t) - -M^V z {^t). 
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and couples to the path x(r) only. With the expression (CIO) for the generating functional and the action (C13), 
we have derived a representation similar to Eq. (Al) with the action ( |A2| ), extended by an additional oscillator in 
z-direction. We identify 



^_L1) — J (^_L2 ^J_l) ) 3x — Jxi jy — Jy 



Thus the auxiliary frequencies Q± (A16) become 



n±(fixi,nj. 2 ) = -\n ±1 ±n ±2 \ 



(C15) 



(C16) 



Inserting the substitutions (C15) into the solution ( A8) for the generating functional in two dimensions and performing 
the usual calculation for a harmonic oscillator with external source |L6l Chaps. 3,5] in z-direction, we obtain an 
intermediate result for the generating functional in three dimensions ( p2[ ): 

^°' X °[j,v] =K*Zn' X °J ^ ex p{^/ l/3rfrV2 ^ r )} dr dr' 3(tr)G^(r,T')3(tr')y 

(C17) 



The partition function follows from Eqs. ( A17 ) and ( A23 ) 



z po,x 0=z po,x 0[0)0] 



npn+/2 h(3Q, {l /2 



smhhpn + /2 smhh(3n_/2 sinhn/3rj||/2 
and G x ° (r, r') is the 3 x 3-matrix of Green functions 

G?i(r,T') G*1(t,t') 



G x °(t,t') = G x S(r,r') G x °(r,r') 







(C18) 



(C19) 







G x o(r,r') 



Except G x °(t, t'), the Green functions are given by the expressions in Eqs. (A21) and (A27) with frequencies (C16). 
The Green function of the pure harmonic oscillator in z-direction 



G x °(t,t') = 



1 



M/3QI V 2 



/h/3Q« cosh ft || ( 



hp/2) 



sinh H/3n»/2 



(C20) 



follows directly from the limit (A22). Since the curr ent J ( C14 ) still depends on time derivatives of V, we have to 
perform some partial integrations in the functional (C17). This is a very extensive but straightforward work and 
thus we only present an instructive example. For that we apply the properties and the time derivatives of the Green 
functions which we presented in Appendix pi Consider the integral 



M 2 



h6 



hp 



(C21) 



occuring in the second exponential of Eq. (C17) with i € {x, y, z}. A partial integration in the r'-integral leads to 



M 2 



2h 2 



K0 



t'=K0 



t'=0 



II 



M 



hp 



I dr I dT'V i ($,T)G-*?(T,T')V i (Z,T'). 



(C22) 



The surface term in the first line vanishes as a consequence of the periodicity of the current ( |C5| ) and the Green 
function (B5). A second partial integration, now in the r-integral, results in 



I=-—s / dr / dT'V^ryG'l^ry^tr') 
in Jo Jo 



2h< 



' J o drj^ dr'V l ^r)'G'::(r,r , )V l (tr')--J^drV^t' 



(C23) 
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Here we have applied the periodicity property of the right-hand derivative of the Green function (B1E), leading to a 
vanishing surface term in this case, too. In the second line, we have used the decomposition (BE) of the double-sided 
differentiated Green fu nction . Note that the last term just cancels the appropriate term in the first expone ntial of the 
right-hand side of Eq. ( C17 ). Eventually, after performing all such partial integrations, we reexpress Eq. (C17) by 



?£°' x ° [j, v] = \£ Zg™ J d 3 £ exp | JLj dT £ dT > T ) H x ° (r, r>) s(£ , r') j 



with six-dimensional sources 



S(C,r) = (j(r),V«,r)) 



(C24) 



(C25) 



and the 6 x 6-matrix H XQ ( r, t ') which has no si gnific ance as long as we have not done the ^-integration. We explicitly 
insert the decomposition (C8) into expression (|C25| ) of t he so urce vector s. Since Vo(£) from Eq. (C7) is constant in 
time, some temporal integrals in the exponential of Eq. ( p24 ) can be calculated and we obtain 



^ X0 [j,v]=A t - h 3 ^< x °exp 



f 1 r hp f hf3 ) 
I2W0 dT Jo dT ' s ( r ) HX °( T ' r ') s ( r ')| 



/ 



cf£ exp 



M 2 ■ M ^ 

2h 2 /3 Ti 2 f3 



hi) 



dr v(r) 



with the new 6-vector 



s(r) = (j(r),v(r)) 



(C26) 



(C27) 



consisting of the original sources j and v only. The Gaussian ^-integral in Eq. (|C26[) can be easily solved and the 
terms appearing from quadratic completion modify the above matrix H xo (r, r'). The final result for the generating 
functional of all position and momentum dependent correlations is given by 



Z%>*° [j, v] = Z^ exp I ^3 ^ dr ^ dr' s(r) (r, r') s(r') \ 



(C28) 



The complete 6 x 6-matrix G Po,x °(r, r') contains all possible Green functions describing position-position, position- 
momentum, and momentum-momentum correlations. As a consequence of separating the fluctuations into those 
perpendicular and parallel to the direction of the magnetic field, all correlations between x, y on the one and z on the 
other hand vanish as well as those for the appropriate momenta. The symmetries for the Green functions and their 
derivatives were investigated in detail in Appendix [b| and lead to a further reduction of the number of significant 
matrix elements. It turns out that only 9 elements are independent of each other. Therefore we can write the matrix 



G X0 ' P0 (r, t') 



/GP£' X0 (t,t') Gpo,*o( T)T /) Gv°?°(t,t>) G™>*°(t,t') 

Gpo,*o (t ' ;T) Gp- X0 (t,t') -G?;f(r,r') Gp- x «(t,t') 

Gp»- X0 (t,t') G p «' X0 (t,t') 

Gg°f(r',r) -GP;f(r',r) ^(r/) G^»(r,r') " 

GPo ; x o(r',r) GP»f(r',r) G^(t',t) G p P^(r,r') 



V 











Gp- X0 (t',t) 











(C29) 



GP- X0 (T,r')y 



The matrix decays into four 3 x 3-blocks, each of the which describing another type of correlation: the upper left 
position-position, the upper right position- momentum (as well as the lower left one), and the lower right momentum- 
momentum correlations. The different elements of the matrix are 



G*r°(T,T' 



= {x{r)x{r'))^ =GH{r,r'), 
= <5(r)y(r')) P J °' X0 =G x »(r,r'), 
^(z(r)z(r'))^' X0 -G x «(r,r'), 



-Am ±1 G x °(r,r'), 



= (i(r)Py(r')) P n °^ = iMG'*°(r, t') + -MSI xi G~%{t, r% 



(C30) 
(C31) 
(C32) 

(C33) 
(C34) 
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(C35) 



G% p I°(t,t>) = (P x (t)P x (t'))%> x ° = -M 2 'G'Z(r,r / )~ l M^ ±1 'G:i(r,r') + ^M^G^ir. r' ) -. (C8.il 



M 
J 



GfX°^ T ') = (PAr)Py(r')) P ^ x " = 1 M 2 'G:Ut,t') - M 2 'G'^(t,t') + -M 2 Q,\ 1 G x °(t,t'), 



Gl^(Ty) = {p z {T)P z {T'm 



M^G'Z^r') 



M 
J' 



(C37) 
(C38) 



where the expectation values are defined by Eq. (3.12). Note that all these Green functions are invariant under time 
translations such that 



GPr o (T,r') = Gp Xo (r-r') 



(C39) 



with fi, v € {x,y, z,p x ,p y ,p z }. 

It is quite instructive to prove that all these Green functions can be decomposed into a quantum statistical and a 
classical part as we did it in Eq. (A21). Since we know that the classical correlation functions do not depend on the 
euclidean time, all derivative terms in Eqs. (C30)-(C38) do not contribute a classical term. We can write each Green 
function 



G p ,1' X0 (t,t') = G^(t,t')-G^ 1 



(C40) 



This relat ion has been alread y checked for Eqs. ( p30| )-( C3^ ) in Appendix |X| The classical contribution is zero in 
Eqs. (C33), (C3J:), and (|C37 ) following from the absence of classical ter ms in derivatives of the Green function s and 
mixed correlations like (A29). It seems surprising that the correlation ( |C34[ ) contains a classical term while (C33) 
possesses none. This is, however, a consequence of the cross product of the orbital angular momentum appearing in 
the action (|C3|) and the explicit classical calculation entails 



G%„ = (x Px ) cl = 0, 



^XPy 



= (*Py> d = 15 



13 n 



2 

±2 



(C41) 



where t he la tter is the subtracted classical term in Eq. (A21) when considering the first two substitutions in (C15). 
In Eq. (C38), the second term is obviously the classical one since 



M 



Gf zPz = (p zPz } ci = j. 



The extraction of the classical terms 



Gi 



M 



= (PxPx) Cl = J [ 1 



n 2 



^±2 ^J_l 



(C42) 



(C43) 



in the case of the Green function G p ° * a (r , t') requires the consideration of the last two terms in Eq. ( |C36 ). Thus we 
have shown that the decomposition (C40) holds for each of the Green functions ( p30| )-( C38D . Note the necessity of 
subtracting the classical terms since they all diverge in the classical limit of high temperatures ((3 — > 0). 
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